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ABSTRACT 



Context. There is mounting evidence for an extra-planar gas layer around the Milky Way disk, similar to the anomalous Hi gas 
detected in a few other galaxies. As much as 10% of the gas may be in this phase. 

Aims. We analyze H i clouds located in the disk-halo interface outside the solar circle to probe the properties of the extra-planar 
H i gas, which is following Galactic rotation. 

Methods. We use the Leiden/Argentine/Bonn (LAB) 21-cm line survey to search for Hi clouds which take part in the rotation of 
the Galactic plane, but are located above the disk layer. Selected regions are mapped with the Effelsberg 100-m telescope. Two of the 
Hi halo clouds are studied in detail for their small scale structure using the Westerbork Synthesis Radio Telescope (WSRT)* and the 
NRAO Very Large Array (VLA)** . 

Results. Data from the 100m telescope allow for the parameterization of 25 distinct Hi halo clouds at Galactocentric radii 10 kpc 
< R <15 kpc and heights 1 kpc< z < 5 kpc. The clouds have a median temperature of 620 K, column densities of Nm~ 10 19 cirT 2 , 
and most of them are surrounded by an extended envelope of warmer H i gas. Interferometer observations for two selected regions 
resolve the Hi clouds into several arc-minute sized cores. These cores show narrow line widths (Awi/2~ 3 kms~'), they have volume 
densities of n > 1.3 cnr 3 , masses up to 24 Mq, and are on average in pressure equilibrium with the surrounding envelopes. Pressures 
and densities fall within the expectations from theoretical phase diagrams (P vs («/,)). The H i cores tend to be unstable if one assumes 
a thermally bistable medium, but are in better agreement with models that predict thermal fragmentation driven by a turbulent flow. 

Key words. Galaxy: halo - radio lines: ISM - ISM: clouds 



1. Introduction 

The Hi gas is a major constituent of the interstellar medium 
(ISM), and it is a well known property of this gas that it settles 
in the Galactic plane. From the very first H i observations it is 
also known that the disk gas is co-rotating with the stellar disk. 
Taking both properties together one may use the gas distribution 
to describe the morphology of the Galactic disk. 

A s yet, there is no sharp boundary for the disk emission. 
lOortl (I1962I) was the first who mentioned this fact. "Well out- 
side the real disk one still finds neutral hydrogen with an average 
density of between 5 and 10 per cent of the intensities one ob- 
serves in the plane". Oort was referring to observations with the 
Dwingeloo telescope, and Shane ( 1967) described this gas later 
as a "galactic envelope", a smooth envelope of neutral hydrogen 
surrounding the spiral structure, following the same Galactic ro- 
tation as the ga s in the plane. F urther di scussion of thi s envelope 
was given by Taka kubol d 1967b and by IShand dl97ll) . but there 
was some concern about a possible contamination by stray radi- 
ation from the antenna diagram of the Dwingeloo telescope. 
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The extra-plana r gas component was also v i sible in the 
IWeaver & Williamsl (119731) survey, and iLockmanl (1 1984 stud- 
ied this feature in some more detail. Supplementing observa- 
tions were made with the NRAO 300-f oot telescope and cam e 
also from the NRAO 140-foot survey bv lBurton & Lisztl dl983l) . 
Lockman argued that his analysis was not affected by stray radi- 
ation. He found for Galactocentric radii 4 ^ R ^ 8 kpc that 13% 
of the H i gas is located outside the disk, extending to z-distances 
of 1 kpc or more and termed this comp onen t an "Hi halo". 

The Bell Labs survey dStark et aljl992l) is only l ittle a ffected 
by instrumental effects, and iLockman & Gehman (1991) used 
this survey to analyze the nature of the vertical Hi gas dis- 
tribution in the direction of the Galactic poles. They proposed 
for the Hi gas a decomposition in several layered structures, 
corresponding to distinct different isothermal cloud populations. 
The scale height for each component results from the pressure 
balance of a cloudy turbulent medium against the gravitational 
potential of the Milky Way. The concept of a layered structure 
of the Hi contains essentially three components: a cold neu- 
tral medium (CNM), a warm neutral medium (W NM) and an 
extra-planar component (Dickey & Lockman 1990) which is of- 
ten called "Lockman Layer". 

The layer concept is based on the average emission from the 
extra-planar H i gas layer which is very faint. The clumpy na- 
ture of the H i gas implies then that extra-planar H i clouds must 
have a low volume filling factor. The layer concept describes 
therefore an ensemble of H i clouds or the probability distribu- 
tion of such objects. First indications for a population of such 
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clumps were found by Simo nsonl d!971l) . These early data came 
from the Dwingeloo telescope. Al most three decades later th e 
Leiden/Dwingeloo survey (LDS. dHartmann & Burtonl 1 19971) ) 
became available and provided a much improved database, more 
sensitive and essentially free of stray radiation. Channel maps 
show num erous clumps and fi laments that are detached from the 
disk, and Kal berla et al.l (1 19981) argued for an extra-planar gas 
layer which can be characterized by a distribution with a veloc- 
ity dispersion of <x = 60 kms" 1 , considerably larg er than the 
dispersion suggested bv lLockman & Gehmanl d 1 99 lb . 

The first high resolution data of the extr a-planar gas layer 
at a beam-width of 9' have been taken by iLockmarj ([2002) 
with the Robert C. Byrd Green Bank Telescope (GBT). These 
observations demonstrated convincingly the nature of the extra- 
planar gas layer as a population of cold clumps with a typical 
mass of 50 Mq. Many of these clumps appear to be surrounded 
by warmer envelop es. A larger sam ple of clouds in the lower 
halo was studied by For d et all (120081) with the Parkes Telescope. 
These clumps are s omewhat l arger and more massive than the 
sample detected by iLockmanl (12002b . This cloud population, lo- 
cated close to R ~ 3.8 kpc, is interpreted as originating from 
a Galactic fountain. Such a model would also explain the high 
kinetic energy which is needed for individual clouds to reach 
large z-distance s. Support for such an interpretation comes from 
IStil et all ( 120061) . They found fast moving clumps in the Galactic 
plane with velocity vectors located within the Galactic plane, 
analogous to fast velocities perpendicular to the plane as sug- 
gested as an explanation for the extra-planar gas layer. 

So far we discussed predominantly Galactocentric distances 
R ;£ 8.5 kpc, where 8.5 kpc is the I.A.U Sun-Galactic cen- 
ter distance, since most of the observations are in this range. 
For a more general description of this phenomenon, in par- 
ticular for the question whether the extra-planar gas layer is 
caused by a founta i n flow , objects at larger distances are needed. 
Kalberla & Dedes (2008) argue that extra-planar gas is present 
even at R ^ 35 kpc. Gas at such distances can hardly originate 
from fountain events. 

Direct evidence for a population of extra -planar H i clouds 
outsid e the Solar circle was first given by Stanimi rovic et al.l 
(2006). Arecibo data in the direction towards the anti-center sug- 
gest that these clouds are not restricted to the inner part of the 
Milky Way disk, which is similar to prehminary resul ts with 
the Effelsberg telescope reported by Kalberl a"et al.l d2005l) . In the 
following we intend to explore the extra-planar gas layer of the 
outer part of the Milky Way in some more detail. Our results are 
based on single-dish observations with the Effelsberg 100-m ra- 
dio telescope and on interferometer observations with the VLA 
and the WSRT array. 

This paper is organized as follows. In Sect.|2]we explain our 
selection criteria for targets that have been mapped with the 100- 
m telescope. Our observations are described in Sect. [3] We dis- 
cuss the properties of the H i cloud sample detected by us with 
the 100-m telescope, and also the results from two targeted inter- 
ferometer observations in Sect. [4] We find evidence for a multi- 
phase structure and compare in Sect. [5] the derived physical pa- 
rameters with theoretical models. Sect.|6]gives our summary and 
conclusion. 

2. Methods 

Extra-planar H i clouds are known to have a patchy distribution. 
Accordingly, a strategy is needed for a successful s earch. If on e 
likes to measure "gas well outside the real disk" (lOortl IT962b . 
the first step is obviously to determine the extension of the disk. 



Next, one needs to search in regions above the disk. For the inner 
part of the Galaxy the approach is easy. The scale height of the 
gas is approximately constant for 3 S= R 5= 8 kpc, the boundary 
between disk and halo is well defined. 13% of the gas resides at 
|z| > 5 00 pc with little fluctuations for 4 ^ R £ 8 kpc (Lockman 
11984 . 

For the outer part of the Milky Way the situation is more 
complex. The gas flares strongly, and in addition the disk is sig- 
nificantly warped. Both cannot be disregarded, and it is neces- 
sary to obtain good estimates for the mid-plane position and 
the scale height of the Hi gas. Two groups have recently in- 
dependently determ ined the shape of t he H i gas distribution, 
iLevine et all d2006l) and lKalberla et all d2007l) . Their results are 
in good agreement, and w e adop t the disk parameters as de- 
rived by Kalberla & Dedes (2008), which were previously also 
used by iKalberla et al.l (120071) for a determination of the aver- 
age extra-planar gas fraction. On average 10% of the Hi gas is 
located outside the disk, the extra-planar gas is well defined for 
8.5 ^ R £ 22 kpc, but tends to incre ase toward small er radii R, 
consistent with the determination by Lockman ( 1984). We con- 
clude that a search for extra-planar H i clouds in the outer part 
of the Milky Way should be promising for R i= 22 kpc if warp 
and flaring are taken into account. 



2.1. Hi halo cloud selection 

We used the lKalberla & Dedesl d2008) model to calculate the ex- 
pected emission t ex (l, b, v) for extra-planar gas and tdiskih b, v) 
respectively for disk gas. The ratio t ex (l, b, v)/tdisk(h b, v) defines 
a probability that a cloud feature, observed at position /, b with 
the velocity v may belong to the extra-planar gas layer. The 
extra-planar gas layer is patchy, and we therefore searched the 
Leiden/Argentine/Bonn (LAB) survey (Kal berla et alj|2005h for 
positions containing weak H i emission that may originate from 
clouds in the lower halo. 

We tested our method in the inner part of the Milky Way and 
recovered those regi ons that have bee n observed previously by 
ILockmanl d2002l) and lFord et alj d2008l) as the ones that are most 
promising for a detection of extra-planar gas clumps. After this 
successful test we applied our search algorithm to R 8.5 kpc. 



2.2. Distance determination 

Quantities directly observable for H i clouds are: column density 
A^hi, angular size of the cloud s and line width Au 1( /2 ■ The cloud 
diameter D, the average spatial volume density (n), pressure P 
and visible mass Mui can be determined only if the distance d 
of the cloud is known. Since the regions probed by us are out- 
side the solar circle, a Milky Way velocity field needs to be used 
to convert the line-of-sight velocity vi sr of a cloud to its dis- 
tance d. We use a m ass model and a rotation curve according to 
IKalberla et al.l d2007l) . which assumes that the halo gas is slightly 
lagging behind the Galactic disk. Assuming co-rotation would 
lead to deviations of 6-22% depending on the region. For the 
iBrand & Blitzl (Il993h rotation curve, differences would amount 
to 15 -30% . Finally , in comparison with the Milky Way model 
from lGom ez (2006) the deviation for the distance determination 
is between -7% and -14%. All together, distances determined by 
us may have typical systematical uncertainties of about 15-20%. 
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Fig. 1. Spectra of Hi clouds detected using the Effelsberg tele- 
scope. Arrows mark the position of the H i clouds in the spec- 
tra. The main galactic line is visible at the right side of the plot. 
The extending wings are seen up to ~-100kms _1 , underlying the 
Hi clouds. 

3. Observations 

3.1. Effelsberg observations 

22 fields, each covering 3°x3°, could be observed with the 100-m 
Effelsberg telescope. All the fields have a longitude of / > 90°, 
and the selection was based on the criteria described in Sect. 
12 Therefore we are confident that the Hi emission is asso- 
ciated with the neutral component of the gaseous halo. The 
observations were done during the period from May 2004 to 
October 2005. We used the AK-90 auto-correlator with two po- 
larizations at a bandwidth of 10 MHz with 2048 channels. This 
results in a channel separation of 4 kHz (velocity separation 
0.84 kms -1 ) and velocity resolution 1.03 km/s. The fields were 
mapped beam-by-beam on a 9' grid. The integration for each 
position was 60 sec. For a system temperature T sys - 21 K this 
implies a sensitivity of cr„=0.1 K. The Effelsberg data were cali- 
brated using the IAU standard position S7 (Kalberla et al. 1982) 
. A first order polynomial was applied to correct the baseline, 
and the st ray radiation contami nation was removed using the 
method bv lKalberla et alj dl980l) . The final result was an image 
cube with a 9' angular resolution and 1 kms -1 velocity reso- 
lution. To verify the observations, the detected H i clouds were 
re-observed using a full sampling (4.5' grid). The configuration 
of the AK-90 auto-correlator was identical. The integration was 
increased to 90sec per position, resulting in a theoretical sensi- 
tivity of cr„ =0.08 K. 

3.2. Synthesis array observations 

After extracting a sample of the Hi halo clouds from the 
Effelsberg data (see Sect. 14. U . follow-up observations were 
made for two of the clouds with the WSRT and the VLA syn- 
thesis arrays. 

The cloud at l,b=116.2°,23.6° at u /s ,.=-68 kms" 1 was ob- 
served with the WSRT array, mapping the region at (J2000) 



a,£=20 h 4 m 22', 82 d 56'36 s in a maxi-short configuration. The in- 
tegration time was 12 hours. A double IF was used with a band- 
width of 2.5 MHz and 1024 channels. This backend configura- 
tion results in a channel separation of 2.5 kHz and a velocity sep- 
aration of 0.5 km s . The data were reduced with the MIRIAE0 
software package. For flux calibration the source 3C286 was 
used, while self-calibration was applied to correct the phase er- 
rors. The dirty cube has a sensitivity of cr„ ~2.25 mJy/beam. 
After applying continuum subtraction, the dirt y cube wa s de- 
convolved using the Clark CLEAN algorithm (Clark1 fT980l) and 
convolved with a Gaussian beam. The final result is an image 
cube with a 60" resolution and 1 kms~' velocity resolution. 

The cloud at l,b=115.0°,+23.9° and u br =-84.50kms- 1 was 
observed with the VLA array mapping the region (J2000) 
(a,6)=20 h 29 m , 82 d 08 ' for six hours in the DnC configu- 
ration. A double IF was used with a bandwidth of 0.78 MHz 
and 256 channels in each IF. With this configuration we have a 
channel separation of 3.05 KHz resulting in a velocity separa- 
tion of 0.64 kms _I . The data were reduced using the NRAO 
Astronomical Image Processing System (AlPSfl The source 
3C286 was used as a flux calibrator and the close-by source 
2344+824 as a phase calibrator. The dirty cube has a sensitiv- 
ity of cr„ ~2.25 mJy/beam. The dirty image wa s de-convolved 
using the Clark CLEAN algorithm dClarkl ll980) after the con- 
tinuum subtraction. The clean components were restored with a 
Gaussian beam, resulting in a clean image with a 60" resolution. 

4. Results 

4.1. Single dish 

The total area covered by the 22 fields observed with the 100-m 
telescope is ~204 deg 2 . Due to the diffuse nature of the clouds 
it was not possible to use automated detection algorithms, so 
the selection was done manually. Application of the criteria dis- 
cussed in Sect. 12.11 yielded 25 objects identified as H i halo 
clouds with line-of-sight velocities v\ sr close to the emission 
of the underlying disk. Fig. Q] shows a few examples of the ob- 
served line spectra. The lines are narrow but superposed on ex- 
tended wings of diffuse Galactic Hi emission. For the 100-m 
telescope, stray radiation effects have been taken into account; 
we are therefore confident that these observed components are 
not caused by instrumental effects. 

Table Q] summarizes our single dish results. From the 
observed cloud position /, b and velocity v/ s , we derive its 
Galactocentric distance R and the height z above the plane 
(Kal berla & Dedesll2008l) . The angular diameter s is the geo- 
metrical mean of the major and the minor axis of the cloud. The 
high brightness temperature Tb of the clouds ensures that the 
effect of the noise is minimal. In some cases where the entry is 
missing in Table Q] the clouds were unresolved and the diam- 
eter could not be constrained due insufficient beam-by-beam 
measurements and because of the blending with the extended 
wings. The spatial diameter is calculated from s making use of 
the known distance. As mentioned in Sec. I2.2l this is one major 
source of uncertainties. Assuming optically thin gas, to obtain 
the line width Avi/2 and the column density Nhi we tried to fit a 
one or two component Gaussian to the average spectrum of the 
cloud. In some clouds such as e.g. 1 16.2+23.6,1 16.7+23.0 this 
was successful. In clouds like e.g. 113.0-12.4 and 113.3-14.1, 
where the extended wing was stronger, this was not possible. 

1 http://www.atnf.csiro.au/computing/software/miriad/ 

2 http://www.aips.nrao.edu/aips_faq.html 
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Table 1. Properties of the observed Hi halo clouds. 
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Note - v lsr is the line of sight velocity in kms -1 , R is the Galactocentric distance, z is the height above mid-plane. N HI is the Hi column density. 
s is the measured angular diameter. D is the spatial diameter. Av 1/2 is the full width half maximum of the line width. T ki „ the kinetic temperature, 
(n) the volume Hi density. P is the pressure of the Hi gas. Mm the visible Hi mass of the cloud. 



In those cases we had to first model the extended wing, then 
subtract it and estimate the Aw 1/2 and the column density of the 
cloud. The observed line width Avi/2 defines the upper limit 
for the kinetic temperature T^ n due to the effect of turbulence. 
For very cold clouds, assuming a spin temperature T s = 80 K, 
we expect the derived temperature not to be strongly biased, 
while for the broader lines the bias is larger. To determine the 
average volume density (n), we assumed a cylindrical shape: 
(n) = NHi/d. Deviations from this shape and the uncertainties 
in the distances determination are the major sources of biases in 
this estimate. The pressure was estimated from P — (n) ■ 7jj n . 
This includes the thermal pressure as well as a turbulence 
component. Finally from the column density Nhi and size of the 
cloud D we estimated its visible mass using M = Nhi * D 2 . 

The Gaussian deconvolution of the clouds showed that most 
of them can best be fitted with a two-component Gaussian. 
This corresponds to a narrow compact component and to a 
broader more extended component. The best two examples 
which are present in this paper are the cloud 116.2+23.6 
(Fig. I3 and 115.0+23.9. All the properties of the narrow 
components are given in Table [TJ except the peak temperature, 
which for cloud 116.2+23.6 is r c =1.8+0.1 K and for the 
cloud 115.0+23.9 is r c =0.73+0.1 K. We estimated a peak 
temperature 7V=0.71+0.1 K, column density AV=25+4- 10 18 
cm -2 and a Afi/2=18+0.5 kms -1 for the broad component of 
the cloud 116.2+23.6. For the cloud 115.0+23.9 the Gaussian 
decomposition gave a peak temperature 7V=0. 8+0.1 K, a col- 
umn density AV=20+4 10 18 cirT 2 and Ayi/ 2 =15+1 kms -1 . Due 
to the presence of the Galactic wings in the emission profiles, 
the size of the broad component is uncertain. Assuming that it 



has a similar dimension as the narrow part, we have a volume 
density of («)=0. 16+0.07 cm -3 and pressure P-\ 151+507 
K ■ cirT 3 for 1 16.2+23.6 and volume density </i>=0. 1 1+0.04 and 
P=535+21 1 K ■ cm -3 for 1 15.0+23.9. 



According to Rohlfs & Wilson (2004), assuming virializa- 
tion, the line widths Avi/2 of the Hi clouds can be used to 
estimate their virial masses. For the Hi clouds in Table[TJthese 
calculations give typical virial masses M V( > ~ 10 4 Mq, which 
are more than two orders of magnitude larger than the visi- 
ble H 1 masses listed there. This comparison indicates that the 
H 1 halo clouds observed with the Effelsberg telescope cannot 
be self-gravitating objects. Therefore an external confinement is 
needed for the clouds not to dispers e . Assu min g the presence 
of a h ot halo according to iPietz et all d!998l) and lKalberla et alj 
(2007), the counterpart could be provided by the envelope and 
the surrounding hot halo medium. But in the cloud 116.2+23.6 
its pressure P=14+8 K-cm" 3 is small in comparison to the en- 
velope pressure of P-l 151+507 K-cm -3 . The same is true for 
the cloud 115.00+23.9, where the pressure for the cloud is 
P=19±7 K-cm -3 in comparison to the pressure of the enve- 
lope of P=535+211 K-cm -3 . In addition, a comparison of the 
Hi cloud pressures P = («) • T with theoretical est imates 
for warm component pressures from Wolfi re et alj d2003l) at the 
given Galactocentric radius range 10 S= R 5= 15 kpc indicate 
a similar trend for the rest of the clouds of the sample. Before 
reaching a firm conclusion that the H 1 halo clouds may be tran- 
sient objects, one can conclude that the clouds are unresolved. 
This implies that estimates for the sizes are upper limits only 
and higher resolution observations are needed. 
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Fig. 2. Left: a) A longitude-velocity Hi brightness temperature map of the cloud 116.20+23.55 taken by the Effelsberg 100-m 
telescope. The cr„ is 0.1K. The color-bar displays the transfer function. The contours are at 0.1, 0.3, 0.6, 0.9, 1, and 2 K. Right:b) 
A latitude-velocity H i brightness temperature map of the cloud 1 16.20+23.55 taken by the Effelsberg 100-m telescope. The cr„ is 
0.1K. The color-bar displays the transfer function. The contours are at the level of 0.1K 0.3K 0.6K 0.9K IK 2K. 



4.2. Synthesis observations 

Two of the H i halo cloud positions from the Effelsberg sam- 
ple were observed at high resolution to get a better constraint 
of the properties of the clouds. The first cloud 116.2+23.6 
was observed with WSRT at a position (J2000) a =20 h 4 m 22', 
82 d 56'36 s and has a line-of-sight velocity u/ !r =-68 kms" 1 , 
which corresponds to a distance of d = 7 kpc (R = 13 kpc, 
Z = 2.5 kpc). It is cold with a line width Ai/i/2=3 kms _I and a 
column density of /Vhi=H ±2 10 18 cm 2 . The second Hi halo 
cloud, 1 15.0+23.9, was observed with VLA, at an offset position 
(J2000) (a,5) =20 h 29 m , 82 d 08' and has a line-of-sight velocity 
vi sr = -84.50 km s _I , which corresponds to a distance of d-10 kpc 
{R = 15 kpc, z — 4 kpc). This cloud shows also a narrow line 
width Aui/2=3 kms -1 and has a column density /Vhi=15 ± 3 
10 18 cm" 2 . 

As seen in Fig. [3^ a collection of compact Hi objects 
was found in the WSRT observations at a i// sr =-68 kms -1 . 
All cores with an angular radius s> 60" are significant at a 
signal-to-noise level of 5 or more. Some unexpected cores were 
found at a line-of-sight velocity of u; jr =-85 kms -1 (Fig[3]b). 
These are associated with an Hi emission that is barely sig- 
nificant in the Effelsberg spectra. Similarly, at the position 
(J2000) (a,<5)=20 h 29 m , 82 d 08', observed with the VLA we find 
a conglomeration of at least 8 cores at a velocity of i/; jr =-84 
kms -1 which are associated with the Hi cloud 115.0+23.9. In 
Fig.|4]a channel map of the VLA observations is given. 

In Table [2] median values for the observed quantities and de- 
rived physical parameters respectively are given for the H i cores 
detected with the WSRT and VLA telescopes. To determine the 
angular sizes and their column densities we produced separated 
0th moment maps, and we fitted a 2D Gaussian. The fitted full 
width half maximum is the angular size s of the cores. From the 
integrated flux we estimated the column density. The line width 
is determined by inspection of the line since it was not possi- 
ble to fit the spectra. Assuming a spherical shape, the average 
volume density is estimated. Finally for the pressure we used a 
similar method as the one described in Sec. 14.11 The Hi cores 
show very narrow lines with a median of Aui/2=3.3 km s _1 and 
At»i/2=4.3 kms -1 for WSRT and VLA respectively, implying 
cold gas. The H i cores are also very compact with a median 
angular size 5=78" for WSRT and 5=88" for VLA, implying 




20:29:23.0 20:28:52.0 20:28:21.0 20:27:51.0 
Right Ascension (J2000) 



Fig. 4. An RA-DEC map of Hi emission of the cloud 
1 15.0+23.9 observed with VLA array centered at a=20 l }29 I ? 1 00 s :, 
( 5=82 d 08'? 1 00 s ., u /jr =-84.50 kms" 1 . Hi emission below the 1- 
cr„ of 0.4K has been blanked out. 

a median linear diameter of D-2.6 pc and D-5 pc respectively. 
For the measured column densities which are given in Table [2] 
the above diameters imply that the cores are a lot denser than 
what was estimated from the Effelsberg observations. For the 
H i cores detected with the WSRT, the median volume density 
is {n)-2. 3+0.8 cirT 3 , while for the Hi cores detected with the 
VLA the median is (n) = 3.6± 1.1cm" 3 . A comparison of the to- 
tal mass observed in the H i cores with the total visible H i mass 
of the "parent" Effelsberg cloud shows that they carry a large 
fraction of the mass. In the case of the WSRT cores they cover 
an area of 15arc-min 2 with a total mass of 33M0, while the VLA 
cores cover an area of 14.55 arcmin 2 and have a total visible 
Hi mass of 107 Mq. The Hi cores have very small filling fac- 
tors of 3%-4%, but carry a significant amount of the mass, up to 
48% in comparison with the Effelsberg observations. The rest of 
the mass is associated most probably with the warmer H i enve- 
lope, which is extended and undetectable for an interferometer. 

The high resolution interferometer observations enable us to 
determine accurately the morphology of the H i clouds observed 
with Effelsberg. They are resolved into individual Hi cores for 
which kinematical pressures could be determined. As seen in 
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Fig. 3. a) A RA-DEC map of Hi emission observed with the WSRT telescope centered at a=2#51™0O s :, 5=83 l ?01™52 s :, u bJ .=-68.72 
kms -1 . Hi emission below the l-cr„ of =0.4K has been blanked out. b) An RA-DEC map of Hi emission of the cloud 1 16.2+23.6 
observed with the WSRT telescope centered at a=20M3™48 s ., 5=83'?04™44 s ., f iir =-84.18 kms -1 . Hi emission below the l-cr„ of 
0.4K has been blanked out. 

Table 2. Derived parameters of the H i cores. 





N H i 




D 


Tkin 


(n) 


P 


M HI 




10 18 cm- 2 


arcmin 


pc 


K 


cirr 3 


K ■ cm -3 


M Q 


Range 


12-32 


58-127 


1.9-5.9 


140-350 


1.3-3.1 


270-980 


1—24 


Median 


16 


78 


2.6 


240 


2.3 


560 


4 


Range 


29-76 


54-124 


2-6 


170-730 


2.2-5.7 


780-2440 


7-32 


Median 


40 


86 


5 


410 


3.6 


1500 


16 



Note - The first two rows correspond to the cloud 116.2+23.6 observed with WSRT. The last two rows correspond to the cloud 115.0+23.9 
observed with VLA. N HI is the observed H i column density of the cores, s is the angular size, D is the spatial diameter. T kin is the kinetic 
temperature, (n) is the volume Hi density. P is the pressure of the Hi gas. M Hl is the visible Hi mass of the core. 



Table [2] the median pressure of P = 560 ± 230 K ern 3 for 
the WSRT sample and the median pressure of P = 1500 ± 490 
K cm 3 for the VLA sample are both comparable within the un- 
certainties with the pressures of the surrounding envelope as esti- 
mated from the Effelsberg observation, which are P-l 151+507 
K-cm 3 , and P=535±211 K-cm~ 3 for the cloud observed with 
WSRT and VLA respectively. 

A comparison of the pressures for cores and envelopes with 
theor etical estimat e s at di fferent Galactocentric radii from Table 
4 of IWolfire et al.l d2003l) shows that the derived quantities fit 
well to a two-phase picture with an approximate pressure equi- 
librium between the cores and the surrounding warm H i en- 
velope. The envelopes as detected in the single dish observa- 
tions (see Fig. |2]l may provide the necessary support to stabi- 
lize the clouds. The envelope with a peak brightness tempera- 
ture of ~0.8 K is too weak to be detectable in our interferome- 
ter data. Additional deep observations would be needed to clar- 
ify the presence and extent of such a component. Nevertheless, 
combining the Effelsberg and WSRT/VLA observations can give 
us an insight into the nature of the ISM in the region of these 
clouds. Taking into account the presence of the very extended 
broad galactic wings, the warm envelope and the small scale 
cold cores which are unresolved with the Effelsberg telescope, 
we reach a sort of hierarchi cal structure of the IS M, similar to 
the turbulent flows found by Audi t & Hennebellel d2005l) . which 
will be discussed later. 

Using equation 47 from iMcKee & Cowid d 19771) we can 
show that the envelope would protect the Hi cores from fast 
evaporation. In the absence of a warm envelope, the mass loss 
for an Hi core embedded in a hot plasma (T~ IQ 6 K) is ~ 
7 • lO^Mg-Myrs- 1 . This implies that a core with an average 



mass of ~ 10 M0 will evaporate in ~ 140Myrs. A core em bed- 
ded in a warm envelope (T ~ 7000 K lWolfire etaf] (120031) ) has 
a mass loss of one magnitude smaller, on the order of ~ 6 • lO -3 
M0-Myrs ', and the cores will evaporate in ~1.7Gyrs. But be- 
cause of internal turbulent motions the clouds will evolve within 
~ 2 • 10 5 yrs. This is a very short time in comparison with the total 
time of their orbit, which is around 100 Myrs, which is found by 
simple ballistic simulation. 



4.3. Comparison between the properties of different samples 
of H 1 clouds 

To better understand the physical properties of the Hi halo 
clouds, but also in order to determine various systematics, it 
is reasonable to compare our clouds with different samples of 
H 1 clouds detected with other telescopes . The samples of clouds 
for a comparison are: a) lLockmanl d2002l) . where a population of 
clouds in the inner Galaxy (R ~ 3.5kpc) was detected with GBT 
The distan ce was determined using t he terminal velocity of the 
sources. b) IStanimirovic et"aT] d2006l) . using the Arecibo 300-m 
radio telescope, detected a number H 1 clouds distinctively sep- 
arated from the Galactic disk towards the anti-center direction. 
Due to high intrinsic uncertainties, kinematic distances were not 
used, and the authors opted to use p ressure equ i librium consid- 
erations to determine the distance. c^ lFord et al.l d2008h . where a 
large number of H 1 clouds were detected within the pilot region 
of the Galactic All-Sky Survey (GASS) with the Parkes 64-m 
telescope. The terminal velocity was used as well to determine 
th e distance , since the clouds are located in the inner Galaxy, 
d) IStil et all d2006l) . where in the VLA Galactic Plane Survey 
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(VGPS) 17 fast moving Hi clouds were detected close to the 
plane which are possibly associated with the halo gas phase. 

An inter-comparison of the main results from these authors is 
given in Table[3] Telescope independent quantities, like Nhi and 
Aui/2, are found to be similar, with an exception of Nhi from 
IStil et alJ d2006h : the higher column densities in this sample can 
be explained by the low altitude of the sample. The other cloud 
properties are similar, which may imply that the different cloud 
samples in the inner and the outer Galaxy may have a similar 
origin. 

The telescope-dependent properties of the samples, s, D, n, 
Mhi, show a different picture. Telescopes with similar resolu- 
tions, like Effelsberg, GBT and Parkes, show similar ranges for 
s, D, (re), Mhi- Taking into account that the source distances 
are fairly similar, it is quite clear that the telescope used in each 
work may introduce a twofold systematic bias: a) Synthesis 
telescopes work as spatial filters, more sensitive to cold compact 
than to warm extended gas. b) The measured angular diameter 
s depends on the convolution of the actual cloud size s t and 
the beam size. As such the s t is overestimated, leading to an 
underestimate of the measured volume density n. For the same 
cloud, telescopes with increased resolution will either resolve 
the cloud or measure s closer to the s t . Therefore the errors in 
the determination of the parameters depending on the angular 
size such as D, n, Mm are smaller, and the measurements are 
closer to the properties of the clouds. This dependence on the 
beam size obviously explains that the Arecibo observations 
(beam resolution ~ 4 ') yield parameters for the clouds closer to 
the ones estimated for the H i cores observed with the WSRT 
and the VLA telescope. 



5. Discussion 

In Sect. 14. 21 we found indications for a pressure equilibrium be- 
tween the Hi cores which constitute the halo clouds and the 
warmer envelopes in which they are p ossibly embedded. I n the 
Galactic disk, it has been shown (see IWolfire et al.ll2003l and 
reference therein) that the two main phases, the cold medium 
phase with T ~ 100K and the warm m edium phase (WNM) with 
T ~ 10 4 K dKulkarni & Heiles| [T987). can co-exist in pressure 
and thermal equilibriu m only in a very nar row range of pres- 
sures and densities. The lWolfire et alJ C2003) model determines 
this range for different Galactocentric distances up to R = 18kpc. 
The model takes into account various observational constraints, 
e.g. dust and metalicities; for a m ore detailed discussion see 
Chap. 2 and 3 of Wo lfire et alJ d2003l) . It is assumed that the main 
heating for the neutral medium originates from the dust grains 
through the FUV of young stars. The cooling of the cold phase 
is mainly due to the fine-structure of the CII line (158/fln), while 
the cooling of the warm phase happens in through the Lya, C II 
(158/vm) and O I (63/zm), with the electron recombination mech- 
anism also playing an important role. It is important to note here 
that the model depends on th e dust-to-gas ratio s, the metalicity 
and the assumed FUV field (IWolfire et al.ll 1995b . 

Since we were able to determine the volume densities of 
the Hi gas (re) and the pressure P for the halo clouds and the 
cores, it is worth co mparing the phase diagrams as estimated by 
IWolfire et al. (2003) with our results. We make the following as- 
sumptions: 

1. In IWolfire et alJ (120031) it is assumed that the total density of 
the hydrogen nucleus isre# = «hi +"m where nm is the spa- 
tially averaged volume Hi density and n#2 is the spatially 



averaged molecular hydrogen density. Up to now no direct 
CO observation of the H i halo clouds has been done, there- 
fore their molecular content is unknown and we assume that 
n H = «m- 

2. The Hi halo clouds detected in this work are located at a 
height of 1 ^ z ^ 5 kpc a bove the plane. According to Eq. 
4 from I Wolfire et all d 19951) . which gives the height depen- 
dence of the FUV field, the FUV for z — 5 kpc is similar to 
the FUV field in the plane (z - 0). Therefore, for dust-to-gas 
ratios and metalicitie s similar to the plan e, we can use the 
phase diagrams from IWolfire et al] (120031) estimated for the 
plane ( z = Okpc). 

In Fig. [5^ we compare the properties of the Hi cores ob- 
served with the VLA/WSRT with the theoretical phase diagrams. 
The envelope estimates in this plot are given from the Effelsberg 
observations described in Sec l4.ll As discussed in Sec. 14.21 the 
H i cores of the VLA have a R ~ 15 kpc, while the WSRT cores 
have an estimated R from 13 kpc up to 15 kpc. As seen in the 
Fig. |5h, at the corresponding R the WSRT cores with v/ sr =- 
68 kms~' are located barely within the range where thermal 
equilibrium is possible and cold gas can exist in a stable phase. 
For the WSRT cores at vi sr — -84km s _I , some are located in the 
thermally unstable region, while others have volume densities 
matching the expectations for the CNM, with smaller pressures 
than the ones expected for the CNM at R=15 kpc. For some of 
the VLA cores in Fig. [5^ we find pressures which are higher than 
expected for cold gas in the equilibrium. 

The derived pressures, plotted in Fig.[5k, may be biased. As 
described in Sec |4.1| they were derived from P — (re) • Tkm- The 
uncertainties are difficult to estimate. The kinetic temperature 
is most probably an upper limit, while the volume density 
(n) is affected by distance uncertainties, beam smearing and 
geometry of the source. Distance uncertainties should be on 
the order of 20%. Beam smearing may cause overestimates 
of the source extension, but more importantly may be biases 
caused by the geo metry of the source. If th ese clouds have a 
sheet-like structure dHeiles & Trol and 2003) we may seriously 
underestimate the volume density and accordingly the pressure. 
Both sources show a very narrow line width. Assuming the 
typical spin temperature in the range of 40 to 80 K, resulting 
biases should be within a factor of a few times. Volume densities 
of the cores may be more seriously biased. Taking this into 
account, the true position of the cores in Fig. [5^ could possibly 
be in better agreement with the cold branch. 

Regarding the position of the cores in the phase diagram, one 
other important detail is that their molecular gas content is un- 
known. Recent absorption experiments report some low column 
density molecular hydrogen in small (O.lpc) and dense clumps 
in the Milky W ay halo with Hi column densities Nh > 10 19 
(iRichterl 12005). Such possible biases imply that the cold gas in 
the Hi cores can be well in thermal equilibrium with the sur- 
rounding warmer gas, detected in the single dish observations, 
which then would play the role of a confining envelope. 

The position of the Hi cores in the phase diagram is ex- 
plained in the above section, using a simplified static thermal 
equilibrium hypothesis between the two components of the ISM. 
Non-thermal effects such as turbulence or magnetic fi elds are not 
taken into account. Recent num erical simulations (Gazo l et al.1 
2005] lAudiL&Hiinebelle 2005) examine the influence of non- 
thermal factors. As a result, the positions of the H i cores in the 
phase diagram (Fig.|5^) and their origins can be explained in the 
context of a dynamical equilibrium under the influence of tur- 
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Table 3. A comparison between the H i halo cloud samples observed with different telescopes. 





Izl 


d 


Nm 


Al»i/2 


s 


D 


(n) 






kpc 


kpc 


ru cm 


km s ' 


i 


pc 


cm ^ 




Effelsberg 


0.9-5.4 


3-11 


0.4-7.4 


2-21.1 


9-27 


8-70 


0.07-1.2 


3-770 


WSRT 






1.2-3.2 


2.5-4.0 


0.9-2.1 


2-6 


1.3-3.1 


1-24 


VLA 






2.9-7.6 


2.8-5.8 


0.9-2.1 


2-6 


2.2-5.7 


7-32 


Lockman (2002) 


0.6-1.2 




0.7-6.3 


5.4-26.3 




19-35 


0.1-0.9 


12-290 


Stanimirovic et al. (2006) 


0.06-0.9 


0.2-3 


1-3 


3-7.6 


6-12 


0.6-8 




0.03-7 


Ford et al. (20081 


0.3-1.7 




0.2-2.2 


5.8-26.2 




46-100 




120-4850 


Stil et al. (2006) 


0.01-0.16 


4.2-7.7 


8-39 


3.4-16.5 


1.8-29.4 


3.4-36.7 


3-24 


9-2500 



Note - The table gives the range of values of the clouds for the following properties: |z| which is the absolute height above the plane (defined in 
the inner part of the disk) in kpc. d the distance in kpc. A'mis the column density of Hi gas in 10"cirT 2 . Avi/2 is the line width in km s" 1 . sis the 
observed angular size in'. D is the diameter in pc. («> is the average volume density in cm 4 . M H i is the visible Hi mass in Mq. 
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Fig. 5. a) Parameters derived from the WSRT & VLA samples. The envelope pressure is determined by the Effelsberg observations. 
The curves in both plots apply to column densities of the order of 10 ■ 10 18 cirT 2 b) A comparison between the Effelsberg sample of 
H i clumps and the phase diagrams depicting thermal pressure P/k vs. hydrogen nucleus density n at different Galactocentric radii 
(IWolfire et al.ll2003h . 



bulence. Audi t & Hennebellel d2005l) examine the influence of 
turbulence in a converging flow of WNM. Under the influence 
of the turbulent velocity field, warm gas is forced out of thermal 
equilibrium into the unstable regime. In this part of the phase 
diagram, the gas forms cold condensations, which are connected 
by less dense structure. 

The more turbulent the gas, the larger the fraction of the 
gas which is driven in the intermediate unstable region of 
the phase diagram ( see fo r comparison in Figs. 3 and 7 of 
lAudit & Hennebellel (12001 1. A very turbulent flow, also part 
of their simulation, generates more complex structures, and in 
addition cold structures are significantly less dense, even at in- 
termediate densities around (n) ~ 5cm -3 . Interestingly, this is 
quite similar to the volume densities we observe in the H i cores 
(Table|2]i. A comparison of the density field and the interferom- 
etry maps (Figs. [3]-@Jl shows a comparable morphology. 

Similar numerical experiments of iGazol et al.l J2005i) study 
the behavior of a bistable gas flow under the influence of a 
turbulent velocity field and model gas that is driven into the 
unstable region. It is shown that as either the effective Mach 
number M or the driving scale increases, a departure of the 
gas from thermal equilibrium is observed, approaching an 
adiabatic behavior. What is interesting in comparison with the 
phase diagrams presented in this work is that a population of 
under-pressured zones is generated in the diffuse gas, while in 
the dense gas over-pressured zones are created. Although, as 



mentioned before, the simulations do not represent an accurate 
model of the ISM, this trend is probably seen in the phase 
diagrams of Fig. [5^. All in all, it seems that the positions of the 
H i cores in the phase diagram agree with the predictions above, 
implying that turbulence strongly affects the state of the halo 
ISM. The Hi halo clouds are probably transient filamentary 
features, cold unstable gas which is continuously condensing 
out of the WNM, a process that is caused by turbulence and 
is stabilized by it. This result fits well to the finding that the 
extra-plana r gas in general is strong ly affected and supported by 
turbulence (Kalb erla & KerpHl998h . 



A comparison of Fig [5^, derived from interferometer 
data, with Fig. |5p, using single dish data only, shows huge 
differences. Parsec sized Hi clumps are unresolved by single 
dish telescopes, estimates for pressures and densities appear 
seriously biased in this case. Interferometers on the other hand 
are insensitive for the extended envelopes which appear to 
surround cold HI cores. The ideal telescope should resolve 
both, the extended features as well as compact cores. Within 
a few years the Australian Square Kilometer Array Pathfinder 
(ASKAP) may be able to satisfy both demands, providing 
a large number of sources in intermediate and high latitudes 
for a comparison with theoretical models of the multi phase ISM. 
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6. Conclusion 

We discussed a population of H i clouds residing in the lower 
halo of the Milky Way, co-rotating with the Galactic disk. 
The sample was observed with the 100-m Effelsberg telescope. 
Search criteria were angular sizes s, the brightness temperatures 
Tb and line width Av\/2 which are considered to be typical for 
halo clumps. The sample includes H i clumps with the following 
properties: 

- they reside in the outer galaxy with Galactocentric radii R 
10 < R < 15 kpc. 

- they belong to the lower halo (0.9 < z < 5.4 kpc). 

- the gas is cold, with a median !).,■„ ~ 600 K and a line width 
Ai»i/2=5.3kms _I . 

- the sample shows a prominent two-component structure. 
Cold H i cores are surrounded by an extended component 
with broad line emission. 

Two of the most prominent H i clouds were observed using 
synthesis arrays, the WSRT and the VLA. These high-resolution 
observations resolve the clouds into a conglomeration of arc- 
minute sized Hi cores. These cores are embedded in a more 
diffuse medium which is detectable only with single dish tele- 
scopes. The cores contain a significant fraction of the Hi mass 
and tend to be in pressure equilibrium with the surrounding en- 
velopes. Taking into account the influence of turbulence onto the 
line widths Avi/2, the median line width values of 3.3 km s'and 
4.3 kms~'observed at the cores implies that the Hi gas is very 
cold. 

Estimating densities and pressures for clumps and surround- 
ing envelopes, we find some scattering but also a reason- 
able agreement with models which predict pressure equilib- 
rium ^nd_j_jnyiW-phase structure caused by thermal instabili- 
ties dWolfire et al.l l2003). The clumps tend to populate unstable 
regions in the phase diagrams, in agreem ent with recent pre- 
dictions of turbulence d riven instabilities dAudit & Henne belle 
l2005HGazol etai1l2005h . 

Compa ring samples ob served with big single dish telescopes 
e.g. GBT (|Lockmanll2002l) . Effelsberg (this work), and Parkes 



T (TLockma 

dFord et al.f l2008), we find similar column densities, peak tem- 
peratures, line widths and masses. Our interferometer observa- 
tions imply that some of the derived parameters may be heavily 
biased if the small scale structure observed by us may be consid- 
ered as typical for H i halo clumps. 
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